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We investigate thermal phase transitions to a valence-bond solid phase in SU(N) Heisenberg 
models with four- or six-body interactions on a square or honeycomb lattice, respectively. In both 
cases, a thermal phase transition occurs that is accompanied by rotational symmetry breaking of 
the lattice. We perform quantum Monte Carlo calculations in order to clarify the critical properties 
of the models. The estimated critical exponents indicate that the universality classes of the square- 
and honeycomb-lattice cases are identical to those of the classical XY model with a Z 4 , symmetry¬ 
breaking field and the 3-state Potts model, respectively. In the square-lattice case, the thermal 
exponent, u, monotonically increases as the system approaches the quantum critical point, while 
the values of the critical exponents, rj and 'f/v, remain constant. From a finite-size scaling analysis, 
we find that the system exhibits weak universality, because the Z 4 symmetry-breaking field is always 
marginal. In contrast, v in the honeycomb-lattice case exhibits a constant value, even in the vicinity 
of the quantum critical point, because the Z 3 field remains relevant in the SU(3) and SU(4) cases. 

PACS numbers: 05.10.Ln, 64.60.De, 75.10.Jm, 75.40.Cx, 75.40.Mg 


I. INTRODUCTION 

The classification of various continuous phase transi¬ 
tions has been successfully discussed from the viewpoint 
of the Landau-Ginzburg-Wilson (LGW) paradigm [l|, Q- 
The essential principles of the paradigm are the clarifica¬ 
tion of (local) order parameters and the characterization 
of breaking symmetries. Recently, the possibility of de- 
confined critical phenomena (DCP) has attracted 
considerable attention as a quantum phase transition 
(QPT) beyond the LGW paradigm. DGP have been pre¬ 
dicted to occur at the QPT point between a magnetically 
ordered phase, such as the Neel phase, and the valence- 
bond solid (VBS) phase in two dimensional (2D) systems. 
Remarkably, this phase transition is continuous, although 
the symmetry group in one phase is not the subset of an¬ 
other phase. The well-known models that are expected to 
exhibit DGP are the generalized Heisenberg models with 
multibody interactions for SU(N) spins namely, SU(N) 
JQm models Gonsiderable effort has been expended 
to numerically determine whether the QPT of this model 
family is of the second order or weak first order, however, 
a satisfactory result has not yet been obtained [mi. 

An interesting aspect of DCP is that the transition may 
occur independently of the lattice geometry ai- In a 
previous study [ll, we evaluated the critical exponent, 
i^QPT, at the QPT point between the Neel and VBS phase 
in SU(N) JQm models on both square and honeycomb 
lattices using quantum Monte Carlo (QMC) calculations. 
From the finite-size scaling (FSS) analysis, we confirmed 
that zzQPT is independent of the lattice geometry but de¬ 
pends on the SU(N) symmetry. This result strongly sug¬ 
gested the presence of DCP in the SU(N) JQm models. 


However, pqpt for the SU(3) models exhibits a system¬ 
atic shift toward the trivial value of pqpt = ^/D{D = 3) 
as the system size increases. Therefore, the possibility of 
a first-order transition remains in the case of SU(3). 

The nature of the QPT point is important in the dis¬ 
cussion of finite-temperature properties, because it can 
strongly affect the topology of the thermal phase diagram 
and also the criticality, as shown in Fig. [T] The SU(N) 
JQm models are expected to exhibit a thermal phase 
transition if the VBS pattern is characterized by sponta¬ 
neous symmetry breaking of the lattice. Thus, consider¬ 
ation of the critical properties of thermal transitions in 
the vicinity of the QPT point may yield a different per¬ 
spective on the possibility of DCP occurring in SU(N) 
JQm models. 

The universality class of the thermal transition has 
been discussed for both SU(2) JQ 2 [l^ and JQ 3 [T^ 
models on the square lattice. The VBS pattern on the 
square lattice is described by a columnar dimer configu¬ 
ration, which is characterized by the spontaneous break¬ 
ing of 7r/2-rotational symmetry around the center of the 
plaquette. Thus, the Z 4 symmetry breaking of the VBS 
order parameter is expected at the critical temperature. 
In the 2D case, several models that exhibit symme¬ 
try breaking exist, such as the Ashkin-Teller model 
including the 4-state Potts model [l^ and the 2D clas¬ 
sical XY spin model with the field (XY-l -^4 model). 
In such models, the critical exponent, 77 , always satisfies 
the condition r\ = 1/4. However, the observed exponent 
77 0.59 of the SU(2) JQ 2 model differs from the ex¬ 

pected value [ 1 ^. In the SU(2) JQ 2 model, the VBS 
order is very weak because the QPT point is located in 
the vicinity of the limit, and the model can only be ex- 
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pressed using the multibody interacting Qm term (the 
dimer limit). To enhance the VBS order, Jin et al. have 
focused on the SU(2) JQ.^ models [l^ . The QMC re¬ 
sults they have obtained [ 1 ^ indicate that the criticality 
is well explained by the Gaussian conformal-field the¬ 
ory with central charge c = 1 ; the thermal exponent, 
v, monotonically increases as the system approaches the 
QPT point, while the following relations between the ex¬ 
ponents, r] = 1/4, 7 / 1 / = 7/4, and j5lv = 1/8, are re¬ 
tained. This is a characteristic aspect of the 2D weak 
Ising universality class [l9j |. and the same behavior has 
also been observed in the 2D XY-I-Z 4 model [20l - [^ . In 
the case of the classical spin model, v monotonically in¬ 
creases as the Z 4 symmetry-breaking field, h^, is sup¬ 
pressed and finally diverges at the XY limit, where the 
Kosterliz-Thouless (KT) transition takes place. The au¬ 
thors in ref. [l^ have observed that an enhancement of 
the U(l) symmetry of the VBS order parameter is ob¬ 
served at close proximity to the transition temperature 
and the QPT point, when the system size is smaller than 
a characteristic length scale. Since it has been noted 
that the emergence of additional U(l) symmetry is an 
important signature of DCP , the numerical result 
in ref. [l^ is consistent with the presence of a deconfined 
critical point in the SU(2) JQ 3 model. However, the ob¬ 
servation of U(l) symmetry in the vicinity of the QPT 
point seems to be natural, because the Z 4 field in the 
classical model is always marginal at a transition tem¬ 
perature and the system becomes the pure XY model 
at the /14 —0 limit [l^. Thus, the emergence of U(l) 
symmetry cannot be regarded as sufficient evidence for 
the presence of a deconfined critical point in this case. 
Since the possibility of a first-order transition has been 
suggested in SU(3) JQ 2 model case , where the same 
Z 4 field is broken, systematic studies of SU(N) symmetry 
are necessary. 

In contrast to the square-lattice case, the nature of the 
symmetry-breaking field is different for the honeycomb- 
lattice case. When the columnar VBS pattern is charac¬ 
terized by 7r/3 rotational symmetry breaking, the corre¬ 
sponding classical model is expected to be the XY-I-Z 3 
model. Since the Z 3 field is relevant in two dimensions, 
the universaliW class is explained by the 2D three-state 
Potts model [2^, and the emergence of the U(l) symme¬ 
try in the VBS order parameter may then be suppressed 
in the vicinity of the QPT. Although this is correct in 
the case of SU(2) spins, the higher SU(N)-symmetric case 
seems to be controversial. The discussion of DCP is based 
on the noncompact complex projective (NCCP^“^) the¬ 
ory with Zk symmetry-breaking fields [3, Q ■ In this the¬ 
ory, although the Z^ symmetry-breaking fields is rele¬ 
vant, it becomes irrelevant as N increases 3 , [ 23 . For 
the SU(2) case, which corresponds to the NCCP^ theory, 
recent QMC results have indicated that the Z 3 field is 
relevant but almost marginal at the QPT [13 . There¬ 
fore, one can expect the first-order transition at the QPT 
point in the SU(2) case and a change of criticality as N in¬ 
creases. This indicates that the criticality of the thermal 


transitions and the topology of the phase diagram are 
determined base on the order of the QPT. If the QPT is 
continuous, as is expected for larger values of N, and the 
system approaches the QPT, whether or not the univer¬ 
sality classes of the thermal transition are affected is a 
nontrivial question. 

Our previous QMC calculations suggest that the same 
criticality exists at the QPT regardless of the lattice ge¬ 
ometry [l3| . This implies that the phase diagram topolo¬ 
gies are identical in both the square- and the honeycomb- 
lattice cases. If one focuses on the most likely and 
simplest case, two scenarios for the thermal phase di¬ 
agram can be expected depending on the order of the 
QPT point: (a) The QPT transition is of the second 
order and the thermal transition is always continuous 
(Fig. (Ha)); and (b) The QPT is a weak first-order tran¬ 
sition and the multicritical point exists at a finite temper¬ 
ature (Fig. [T](b)). When scenario (b) occurs, we expect 
to observe crossover behavior and for v to change to the 
trivial value, v = 1/D{D = 2). From the above dis¬ 
cussion, the importance of calculating the thermal phase 
diagram for different values of N and various lattice ge¬ 
ometries with high accuracy is apparent. Further, such 
calculations can allow us to consider the possibility of 
the DCP scenario in the SU(N) JQm models. Thus, in 
this paper, we systematically study the thermal phase 
transitions of the JQ2 model on the square lattice and 
the JQa model on the honeycomb lattice for SU(3) and 
SU(4) spins. 

The layout of this paper is as follows. In Sec. II, we 
study the thermal transition of the SU(N) JQm model. 
We begin by introducing the model details and the order 
parameters evaluated in the QMC computations. In Sec. 
Ill, we present the results of the finite-size scaling analy¬ 
sis for the obtained numerical data. The criticality of the 
thermal transition is discussed for the square-lattice and 
the honeycomb-lattice cases. Then, we discuss possible 
scenarios for the QPT of both models from the perspec¬ 
tive of the thermal phase diagram. Finally, we summarize 
our results in Sec. IV. 


II. MODELS AND METHOD 

We consider the SU(N) JQ 2 model on the square lat¬ 
tice and the SU(N) JQ 3 model on the honeycomb lat¬ 
tice. Both models are simply expressed by the color- 
singlet-projection operator, , which is defined as P^ = 

■«"here is the SU(N) spin 

generator and is its conjugate. The model Hamilto¬ 
nian can be expressed as 

n = -jY,Pj-Q2 Y. ( 1 ) 

(ij) {i3){kl) 

for the square-lattice case and 

n=-jYPn-Qi Y P^PhlPmn, ( 2 ) 

(b) {i 3 ){kl){mn) 
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SU(N) JQi JQs 

2 Ac = 0.042 Ac = 0.456 

3 Ac = 0.665 Ac = 0.796 

4 Ac = 0.917 Ac = 0.985 


TABLE I: Critical points of SU(N) JQm models. Ac is the 
critical value of the coupling ratio defined as A = J/( J + Qm), 
where J and Qm are the coupling constants. All values are 
given in ref. [Q. 


FIG. 1: Schematic phase diagram and renormalization flow. 
The thick solid (dashed) curves correspond to the second 
(first) order transition. The horizontal axis, A, is the coupling 
ratio of the Heisenberg term, J, and the multibody interac¬ 
tion term Qm- The open square represents a discontinuons 
transition. Each solid circle denotes a fixed point, such as 
the 2D Ising, three-state Potts, and multicritical fixed points. 
The coordination origin corresponds to the low-temperature 
fixed point. All arrows indicate renormalization flows, (a) 
DCP scenario and (b) first-order transition scenario. 


for the honeycomb-lattice case, where (ij) indicates the 



FIG. 2: (Color online) (a) Golor-singlet projection operator 
on a bond. The bold ellipsoids denote a color-singlet dimer 
state and correspond to PijS. (b) Projection operators for Q 2 
and Qs terms, (c) Coordination index, p. 


ponent of the projection operator, z is the coordination 
number of a lattice, and represents the neighboring 
site of r in the fi direction, respectively (see Fig. [5] (b)). 
From lE'.r, the VBS order parameter, which is defined as 
'll = After was evaluated, we obtained 

further quantities: the Binder ratio Br = 
the VBS correlation function, C(r) = the cor¬ 
relation ratio Cr = the correlation length, 

^ ~ static structure fac¬ 
tor, S{Q) = exp[-IQ(r - r')](5'.rT'r.). Here, 

AQ denotes the distance between the order wave-vector, 
Qc — 0, and the nearest-neighbor positions, (0,27r/Ly) 
or (27r/La;,0). 

In this paper, we discuss the thermal transition criti¬ 
cality by changing the coupling constants, J and Qm- It 
is convenient to introduce a length scale associated with 
the distance from the QPT point, where the ground state 
changes from the Neel state to the VBS state. The QPT 
points were previously evaluated in ref. [l^ and are sum¬ 
marized in tablelD The coupling ratio, X = J/{J + Qm), 
of the QPT point depends strongly on the lattice geome¬ 
try, and also on the degree of freedom of the SU(N) spin. 
Therefore, we introduce a normalized coupling constant 
that is defined as A = A/Ac, where Ac is the critical value 
at the QPT point. From this definition, one can easily 
see that A = 0 and I correspond to the dimer limit and 
the QPT point, respectively. 


III. NUMERICAL RESULTS AND FINITE-SIZE 
SCALING ANALYSIS 


nearest-neighbor sites. The summation for the Qm terms 
runs over all pairs without breaking the rotational sym¬ 
metry of the lattice, as illustrated in Fig. [2] Since the 
present lattices are bipartite, the fundamental (conju¬ 
gate) representation is adapted for the SU(N) spins on 
A(B) sites. 

For the Hamiltonians ([T]) and ([2]), we performed QMC 
calculations up to L = 256 for the square- and L = 132 
for the honeycomb-lattice cases, respectively. (The num¬ 
ber of sites, Af, corresponds to Af = L'^ and Af = 2L^, 
respectively.) The QMC code used here is based on the 
massively parallelized Loop algorithm provided in 
the ALPS project code [^. In the computations, we 
measured the VBS amplitude, which is defined as = 
exp[^^].Pr,r^, where Pr,r^ is the diagonal com- 


In Figs. [3]and|31 we show the temperature dependence 
of Cr, Br, and S{Qc) at A = 0.5, which is the middle 
distance between the QPT point and the dimer limit. 
Since clear crosses are always observed for 0 < A < 1 as 
the temperature decreases, the thermal transition from 
the paramagnetic to the VBS phase is expected to be of 
the second order. 

To discuss the universality class, we performed a FSS 
analysis of f, Cr, Br, and S(Qc), assuming the scaling 
forms, ^/L ~ g^[Ly^{T - T,)], Cr ^ gcALy^iT - T,)], 
Br gsAL^^iT - TA, and S(Qc)L-i - - 

Tc)], where yt = and gx[x] is a scaling function. 
We applied the Bayesian scaling analysis to the FSS 
analysis of larger system size sets and estimated the val- 






















4 






CQ 


JJJ* 


hs 

to 


0.50 

0 

2 

1.50 


1.00 

0.50 

0 

0.08 

0.06 

0.04 

0.02 

0 






0.2 


0.3 


ksTI(J+Q2) 


FIG. 3: (Color online) Temperature dependence of Cr^ Br, 
^/L, and S{Qc)L~~ in the SU(3) square-lattice model at A = 
0.5. 
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FIG. 4: (Color online) Temperature dependence of Cr, Br, 
f/L, and S{Qc)L~~ for the SU(4) honeycomb-lattice model 
at A = 0.5. 


ues as follows. First, the critical temperature, T^, and 
yt were evaluated from ^ and Cr (or ^ and the Binder 
ratio Br), because their scaling forms contain only two 
variables, Tc and yt- Both Tc and yt were optimized si¬ 
multaneously from the ^ and Cr data set. In detail, we 
evaluated and yt for several data sets labeled Lmax 
that include four different system sizes, for example. 




FIG. 5: (Color online) (a) Critical temperature and (b) 
renormalization group eigenvalue, yt, for temperature in the 
square-lattice case. The open squares (circles) are the SU(3) 
(SU(4)) results, yt is estimated by extrapolation to the ther¬ 
modynamic limit, A = 0 corresponds to the dimer limit, where 
J — 0, and A = 1 is the QPT point. 


Lmax = 128 includes L = {48,64,96,128}, Lmax = 96 
includes L = {48,60,72,96}, and so on. Since appar¬ 
ent system-size dependence is observed for A > 0, we 
evaluated the extrapolated values of Tc and yt in the 
limit Lniax oo from the large-system sets. (One exam¬ 
ple of this size dependence is the result at A = 0.15 for 
SU(3) shown in Fig. [T]) After we obtained Tc and yt for 
the thermodynamic limit, rj and y/iz were independently 
from the correlation function, C{r), and S{Qc)- 

We summarize the estimated yt{= and Tc in Fig. 
[^for the square-lattice case (the SU(N) JQ^ model). In 
the square-lattice model, yt {v) monotonically decreases 
(increases) as the system approaches the quantum critical 
point, in both the SU(3) and the SU(4) cases. In contrast 
to yt, we observe that rj and y/jz take constant values for 
A < 0.97. Figures [6] (a) and (b) show the A dependence 
of the effective y estimated from the assumption that 
C{R) ~ R~'^. From Figs. [S] (a) and (b), it is apparent 
that y clearly crosses y = 1/4 at critical temperatures 
within the error bars. In the same manner, the effective 
is estimated from the form, S{Qc = 0) ^ L^/'^. Fig¬ 
ures El (c) and (d) present y/jz evaluated from the data 
for L > 96. We can confirm from Fig. IHlthat ^jv crosses 
the value 7/4 at critical temperatures. Thus we conclude 
that y and y/zz satisfy 77 = 1/4 and y/z/ = 7/4 at critical 
temperatures, within the error bars. The obtained expo¬ 
nents are the same as those of the 2D Ising universality 
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FIG. 6 : (Color online) A dependence of effective r; and 7 / 2 / 
of SU(N) JQ 2 models. All values were evaluated from the 
assumptions, C{R — Z//%/2)|t~Tc ~ L~'^ and S{Qc)\t^Tc ~ 

_7 

L , which are approximately satisfied in the vicinity of the 
critical temperatures. The vertical colored lines are critical 
temperatures and the black horizontal lines correspond to the 
values of the exponents for the 2D Ising universality class 
(7 = 1/4 and 7 / 1 / = 7/4)- 
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FIG. 7: (Color online) System-size dependence of yt estimated 
from Br. All yt values were evaluated from the Bayesian scal¬ 
ing analysis for several data sets labeled Lmax (see text). The 
extrapolated values are expected to be those in the thermo¬ 
dynamic limit. The dotted line is a guide for the eye. The 
inset is the same result plotted as a function of /f/Tmax. 


class, yt (v) itself varies depending on A, but the other 
exponents, such as 77 and j/v, are constant. This behav¬ 
ior is known as the 2 D Ising weak universali ty Il9ll and 
is consistent with the results reported in ref. |16l | for the 
SU(2) JQ 3 model. 

To approach the QPT point from the finite tempera¬ 
ture region, we performed these calculations at very low 
fixed temperatures by varying A. With limited system 
size, we observed an apparent increase in yt. However, 
as we discuss below, this is due to crossover from the 
mean-field type behavior to the true asymptotic behav¬ 
ior, and should not be taken as an evidence suggest¬ 
ing a first-order transition. (This is a slightly confusing 
point since the mean-field value for yt = 2 happens to 
be equal to the expected value for the first-order transi¬ 
tion in two dimensions.) Figure [7] shows the system size 
dependence of yt at ksTc/i J -\- Q 2 ) = 1 / 20 , 1/32, 1/64, 


and 1/256 for the SU(3) case. Each value of yt is the 
FSS result of Br for several data sets labeled Tmax that 
include three different system sizes; for example, in the 
square-lattice case, Tmax = 48 includes L — {24,32,48}, 
imax = 72 contains L = (48,64,72} and so on. At 
ksTc/i J + Q 2 ) = 1/20, yt systematically decreases as 
Lmax increases and takes the approximate value yt ~ 0.23 
in the thermodynamic limit. However, in the lower tem¬ 
perature region, we observe that yt exhibits a crossover 
from the mean-field value; the data for small Tmax indi¬ 
cate yt ~ 2(= i), but yt decreases suddenly when the 
system size becomes larger than a characteristic length, 
Lc- In the case of ksTc/iJ + Q 2 ) = 1/32 and 1/64, we 
estimated Lc ^ 72 and Lc ^ 192, respectively. However, 
when kBTc/{J + < 52 ) = 1/256, we obtained a data with 
the exponents of the mean-field value in both the SU(3) 
and SU(4) case. Therefore, we can obtain the correct 
values from the data for L > Lc va. the FSS analysis, 
while we estimate the mean-field values from the L < Lc 
data. This Lc is natively related to the development 
of the correlation length along the imaginary-time direc¬ 
tion, ^ 2 -; the thermal criticality can be observed after 
approximately exceeds the inverse temperature, /3. (It 
is expected that Lc ^ where a is an unknown 

constant.) In the present case, the correlation along the 
real space direction is well developed for < P- Thus, 
the system can be described by an effective model with 
long-range interactions. Similar crossover is observed for 
the critical exponent yt (= in the 2D Ising models 
with long-range interactions [28l |. In the Ising model, yt 
depends on the ratio between the interaction range and 
system size. When the interacting range is significantly 
larger than the system size, mean-field-type behavior is 
observed. From the extrapolated results for yt, it can be 
stated that the universality class of the thermal transi¬ 
tion for Tc > 1/64 is explained by that of the 2D classical 
XY-I-Z 4 model, and is therefore the weak 2D Ising uni¬ 
versality class. 

Next, we focus on the criticality of the SU(N) JQs 
model on the honeycomb lattice. In Fig. [51 we summa¬ 
rize yt=i'~^ and Tc for the SU(3) and SU(4) JQ 3 model. 
We find that yt=()/5 {v = 5/6) is well satisfied even in 
the vicinity of the QPT limit of A = 1 and that the size 
dependence of yt is quite small for A < 0.95. This value 
is consistent with that of the 2D three-state Potts univer¬ 
sality. In Fig. in we show 77 and 7/72 that were estimated 
in the same manner as in the square-lattice case. The 2D 
three-state Potts universality is also confirmed directly; 
77 = 4/15 and 7/72 = 26/15 are satisfied at the critical 
temperatures within error bars. The present columnar 
VBS pattern is characterized by the 7 r/ 3 -rotational sym¬ 
metry breaking, reflecting the honeycomb-lattice back¬ 
ground. Thus, it is expected that the related classical 
model with the same universality class is the 2D XY-l -^3 
model. Since the held is strongly relevant in two di¬ 
mensions, the exponents are not affected by the coupling 
constants J and Q^. 

The fact that the Z^, held is relevant may help us dis- 
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FIG. 8 : (Color online) (a) Critical temperature and (b) 
renormalization group eigenvalue for temperature in the 
honeycomb-lattice case, yt is also evaluated from extrapo¬ 
lation to the thermodynamic limit. The inset of (b) is the 
system-size dependence of yt for the results obtained from 
the fixed-temperature calculations at kBTc/{J -|- Q3) = 1/20 
and 1/64 (see in text). The label tc in the figure means 
ksTc/iJ + Qs). The open squares (circles) are the results 
for the SU(3) (SU(4)) spins. The black dotted line is the 
value of the 2D three-state Potts case, yt = 6/5. 


CUSS the possibility of DCP occurring. If the present 
honeycomb JQ 3 model can be well mapped onto the 2D 
three-state Potts model and the change in the coupling 
ratio can be regarded as the variation of certain parame¬ 
ters, for example, the transverse field in the conventional 
2D Ising model, the criticality in the QPT limit is ex¬ 
plained by the 3D three-state Potts model. In that case, 
the QPT should exhibit a weak first-order transition m- 
and the first-order transition line should extend in the 
finite-temperature region. The length of the first-order 
transition line may be finite but is too short to be ob¬ 
served (see Fig. [D (b)). This means that the value of iy 
should change from the 2D three-state Potts value to the 
trivial value o^ ly — 1/D {D = 2) via the value at the 
multicritical fixed point. However, such crossover be¬ 
havior is not observed when the system approaches the 
QPT point. We also performed the fixed-temperature 
calculations for the honeycomb-lattice case. When we 
vary A for fixed ksT/^J -h Qs) = 1/64, where the crit¬ 
ical point corresponds to A ^ 0.99, the 2D three-state 
Potts universality is still observed. Figure flOl shows the 
FSS results for the SU(3) case. We obtain data collapse 
for L > 80 if we set the critical exponents to those of 
the 2D three-state Potts universality. In the case of the 
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FIG. 9: (Color online) A dependence of estimated 77 and 7 / 7 / 
for JQz model on honeycomb lattice, (a) and (b) ((c) and (d)) 
are rj ( 7 /i^) results for the SU(3) and SU(4) cases, respectively. 
The vertical colored lines denote critical temperatures and the 
black horizontal lines are critical exponents for the 2D three- 
state Potts model, rj = 4/15 and 'y/u = 26/15. 
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FIG. 10: (Color online) Finite-size scaling analysis for SU(3) 
honeycomb-lattice case at fcsTc/(J + Q 3 ) = 1/64. (a) Corre¬ 
lation length and (b) static structure factor. 


honeycomb-lattice model, we expect that the crossover 
behavior from the mean-field theory exists for L < 80, 
but it is very weak. Therefore, it is difficult to iden¬ 
tify the conventional system-size dependence. This re¬ 
sult indicates that the development of Lc is relatively 
suppressed in comparison with the square-lattice case at 
the same temperature. 

The obtained thermal phase diagram for A < 0.99 sup¬ 
ports the possibility of scenario (a) in Fig. [U because 
it seems unlikely that v will approach the trivial value 
of 1/D in both the square-lattice and the honeycomb- 
lattice cases. If the scenario (b) occurs, the multicriti¬ 
cal point should exist at quite a low temperature, i.e., 
A:bT/(J -I- Qm) < 0(10“^). This is still consistent with 
our previous discussion of the QPT point 0; a system¬ 
atic increase in izqpt towards the trivial value is observed 
for L > 128 in the SU(3) square-lattice model. If the dy¬ 
namical exponent for the DCP is unity, ksT/[J-\-Qm) < 
0(10“^) corresponds to the length scale, L > O(IO^). 
This implies that the correlation length is very large and 
almost diverging. 
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IV. SUMMARY 

In this paper, we have investigated the thermal transi¬ 
tions of JQ 2 models on the square lattice and JQ 3 mod¬ 
els on the honeycomb lattice for SU(3) and SU(4) spins. 
We have found that the criticality of the SU(N) square- 
lattice model is well explained by the 2D weak Ising uni¬ 
versality class in both the SU(3) and SU(4) cases, which 
is in agreement with Jin and Sandvik’s result [l6l| for 
the SU(2) JQ3 model. The thermal exponent, v, mono- 
tonically increases as the system approaches the QPT 
limit, and the decrease in v that should occur if 1 / even¬ 
tually reaches its first-order transition value of 1/D has 
not been observed. Thus, the first-order transition ap¬ 
pears to be less likely for ksTc/iJ + Qm) > 0(10“^). In 
the honeycomb-lattice case, reflecting the fact that the 
Z 3 field is strongly relevant, i' always exhibits the 2D 
three-state Potts value. From the obtained results, we 
have discussed possible scenarios for the thermal phase 
diagram. If the first-order transition occurs, we may ob¬ 
serve critical behaviors with strong system-size correc¬ 
tions. However for kBTc/{J + Qm) > 1/64, cross-over 


behavior is not observed clearly in our results. To de¬ 
termine the thermal phase diagram (a) or (b) occurring 
in the present models, the numerical calculations for ex¬ 
tremely large system sizes are required, because the dras¬ 
tic development of Lc is expected in the vicinity of the 
QPT. 
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